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ABSTRACT 

We present the first successful simulation of a neutrino-driven supernova explosion in three dimensions 
(3D), using the Prometheus-Vertex code with an axis-free Yin-Yang grid and a sophisticated treatment of 
three-fiavor, energy-dependent neutrino transport. The progenitor is a nonrotating, zero-metallicity 9.6 Mq star 
with an iron core. While in spherical symmetry outward shock acceleration sets in later than 300 ms after 
bounce, a successful explosion starts at ~ 130 ms postbounce in two dimensions (2D). The 3D model explodes 
at about the same time but with faster shock expansion than in 2D and a more quickly increasing and roughly 
10% higher explosion energy of >10^^ erg. The more favorable explosion conditions in 3D are explained 
by lower temperatures and thus reduced neutrino emission in the cooling layer below the gain radius. This 
moves the gain radius inward and leads to a bigger mass in the gain layer, whose larger recombination energy 
boosts the explosion energy in 3D. These differences are caused by less coherent, less massive, and less rapid 
convective downdrafts associated with postshock convection in 3D. The less violent impact of these accretion 
downflows in the cooling layer produces less shock heating and therefore diminishes energy losses by neutrino 
emission. We thus have, for the first time, identified a reduced mass accretion rate, lower infall velocities, 
and a smaller surface filling factor of convective downdrafts as consequences of 3D postshock turbulence that 
facilitate neutrino-driven explosions and strengthen them compared to the 2D case. 

Subject headings: supemovae: general — hydrodynamics — instabilities — neutrinos 


1. INTRODUCTION 

Modeling the core-collapse supernova (SN) mechanism in 
three dimensions (3D) is still in its infancy. Generaliza¬ 
tion from spherical symmetry (ID) to axial symmetry (2D) 
introduces nonradial flows and hydrodynamic instabilities 
like conv ection and the standi ng accretion shock instability 
(“SASI”; [Blondin et al.||2003| ) in the neutrino-heated post¬ 
shock layer, which have been recognized as helpful for the 
explosion due to improved neutrino-heating conditi ons and 
buoyancy or turbulent pressure behind the shock (e.g.,|Herant 


plosions in 3D, because it destroys the biggest convective 


et ai. 1 199^ [Burrows et al.|l99511 Janka & Muller|1996tlM^ 

>hy & Burrows||2008t [Murphy et al.||2013t |Marek & Ja^ 


009[|Muller et al.||^01^a|b||^bl j[|Mfiller & Janka||^0l4a|b[ 

Suwa et al.||2Q131 IBruenn et al.||2()131 |2Q141 [Takiwaki et al.| 


2014[[Couch & Qtt|2014| ). 

Explosions of 2D models, however, start late and their en¬ 


ergies tend t o be fairly low (except for those of [Bruenn et al. 
2013 2014| whose results still need better explanation and 


confirmation by detailed comparisons with other recent sim¬ 
ulations). 3D effects were hoped to improve the situation, 
but a first optimistic report ba sed on parametrically t riggered 
neutrino-driven explosions by jNordhaus et al. (|2010|) was not 
supported by subsequent works (jHanke et al.|2012 Couchj 


2013t [Couch & Q’ConnorPOT^jTakiwald et ai.|2012 |2014t 


Mezzacappa et al. 

12015). 

Dolence et al. ( 

2013J and Burrows 

et al.[ ([2012[) confessed an error in [Norc 

lhaus et al.[ ([2010[), 


nevertheless they still claimed a remaining positive effect in 
3D. Turbulent fragmentation by energy cascading from large 
to small scales, however, seems to disfavor and delay ex- 
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plumes, which were 

recognized as helful for 2D explosions 

dHanke et al.||20121 

Couch||2013t|Couch & O’Connor|2014| 

[Dolence et al.[[2013 

). Therefore [Abdikamalov et al.[}2014|) 


are worried about even less favorable explosion conditions 
when resolution s hortcomings of all current 3 D models can 
be overcome, and [Couch & Qttj ( [2013[ [2014j) advocate pre¬ 


collapse progenitor-core asymmetries as possible solution of 
the dilemma. 

Here we present the first successful 3D simulation of a 
neutrino-driven SN explosion of a 9.6 M© iron-core star com¬ 
puted fully self-consistently with the neutrino-hydrodynamics 
code Prometheus-Vertex. We identify, for the first time, con¬ 
sequences of 3D turbulence in the convective gain layer that 
enhance the explosion energy and accelerate shock expansion 
in 3D relative to 2D. While so far in more massive progenitors 
3D turbulence appears less supportive for the ini tiation of ex¬ 


plosions than 2D flows (cf. r eferences above and [Hanke et al.[ 
[2013t [Thmborra et ar[[2014[ for 11.2,20,27 M© models), we 
investigate here whether 3D can strengthen the explosion af¬ 
ter the onset of shock runaway. The considered 9.6 M© model 
offers an ideal case for this study because it explodes at the 
same time in 2D and 3D. We briefly describe our numerical 
approach in Sect. discuss our results in Sect. and con¬ 
clude in Sect.m 

2. NUMERICAL SETUP AND PROGENITOR MODEL 

We performed ID, 2D, and full (4;r) 3D simulations of 
a nonrotating, zero-metallicity 9.6 M©, iron-cor e progenitor 
provided by A. Heger (private communication; [Woosley & 


Heger[2015b and previously investigated in 2D by [Muller et al. 

iloiil - 

We used the Prometheus-Vertex hydrodynamics code with 
three-fiavor, energy-dependent, ray-by-ray-plus (RbR-f) neu¬ 
trino transport including the full set of neutrino reactions and 
microphysics ( Rampp & Janka[|2002[ [Buras et al.[[2006 ) ap- 
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Isosurfaces: Entropy/Nucleon Isosurfaces: Entropy/Nucleon Isosurfaces: Entropy/Nucleon 

Colors: Radial Velocity, le9 cm/s Colors: Radial Velocity, le9 cm/s Colors: Radial Velocity, le9 cm/s 
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Fig. 1.— 3D iso-entropy surfaces at 90, 170, and 350 ms after bounce. Colors represent radial velocities. The supernova shock is visible by a thin surrounding 
line, the proto-neutron star by a whitish iso-density surface of 10^^ g cm“^. The yardstick indicates the rapidly growing volume. 


plied i n 3D also by |Hanke et al. ( 2Q13| ) and |Tamborra et aL 
( 2014), in particular the high-density equation of state of |Lat- 
|timer & Swesty| (|19911) with a nuclear incompressibility of 
K = 220 MeV. Our simulations were conducted with a ID 
gravity potential (which is unproblematic for nearly sp herical 
explosions) including general relativistic corrections ( |Marek| 
|et al. 112006] ). For the first time we employed the newly im ¬ 
plemented axis-free Yin-Yang grid ([Kageyama & Sato|2004} . 
The implementation followed [Wongwathanarat et al.| ( |2010] ) 
and posed no particular problems for the RbR+ transport. 
Conservation laws are globally fulfilled with an accuracy of 
-10“^ over several 100 ms for our angular resolution of 2°. 
The radial grid had a reflecting boundary condition at the co¬ 
ordinate center and an infiow condition at the outer boundary 
of 10^ cm. It had 400 nonequidistant zones initially and was 
refined in steps up to >600 zones, providing an increasingly 
better resolution of Ar/r ~ 0.01... 0.004 around the gain ra¬ 
dius. For the neutrino transport 12 geometrically spaced en¬ 
ergy bins with an upper bound of 380 MeV were used. 


3. SIMULATION RESULTS 

The relatively steep density gradient above the iron core 
enables a neutrino-driven explosion of the 9.6 M© star at 
roughly 300 ms after bounce even in ID (left panels in Fig.[^, 
similar to (but more delayed than) oxygen-neon-magnesi um- 
core progenitors exploding as “electron-capture SNe’ ’ ( |Ki-| 
[taura et ^.|2006|panka et al.|2008|[Fischer et al.|2010| ). The 

energy of such a late explosion, however, remains low, only 
£expi ~ 2 X 10"^^ erg (Fig. S- because it is provided mainly 
by the recombination of free nucl eons to a-particles a nd iron- 
group nuclei in the gain layer (see |Scheck et al.|2006| ) and the 
corresponding mass at fpb > 300 ms is only 2.2x10“^ M© in 
ID. Of course, the subsequent neutrino-driven wind from the 
proto-neutron star (PNS) will increase the power of the ex¬ 
plosion. £’expi denotes the instantaneous “diagnostic energy”, 
which is defined by 


^expl 


-L 


dV petot ■ 


( 1 ) 


etot>0,postshock 


Here, p is the density, and the volume integration is performed 
over the postshock region where the total specific energy, 

^tot — ^ + O -f j^^bind( F^) “ ^bindj ? (2) 

is positive, with e, ^|vp, and O being the specific internal, 
kinetic, and (Newtonian) gravitational energies. The brack¬ 
eted term expresses the difference between the nuclear bind¬ 
ing energies per unit mass of all nucleons finally recombined 
to iron-group nuclei and for the nuclear composition at a 
given time. It therefore accounts for the maximum release 
of nuclear binding energy and corresponds to an upper limit 
of £’exp, while omitting this term yields a lower bound on 
£’exp (red and blue lines, respectively, in bottom left panel of 
Fig.^. The binding energy of the stellar layers ahead of the 
shocK plays only a minor role in the energy budget of the ex¬ 
plosion because at ^pb = 400 ms it is only -3.5 x 10"^^ erg in 
the ID simulation and even only -9.5 x 10"^^ erg in our multi¬ 
dimensional simulations. 

3.1. Explosion Dynamics and Properties: 2D 3D 

In 3D convective overturn in the neutrino-heated post¬ 
shock layer develops at ^pb > 70 ms, showing the well-known 
Rayleigh-Taylor mushrooms (Fi g. | T| left) and increasing non- 
radial velocities (colors in Fig. The postshock convec¬ 
tion reaches its maximum activity between about 100 ms and 
200 ms after bounce with a prominent maximum at ~ 120 ms, 
at which time the shock starts its accelerated expansion (cf. 
Figs. and ^ and the postshock flow becomes highly turbu¬ 
lent (Fig. [^middle). At ^pb >200 ms the shock and post¬ 
shock matter expand with -25,000 km s“^ and the ejecta at¬ 
tain a nearly self-similar structure (Fig. right). This is the 
time when convection and turbulence behind the shock be¬ 
come weaker again (less intense red in Fig.|^ and mass shells 
leaving the PNS surface indicate that the compact remnant be¬ 
gins to lose material in the low-density, high-entropy baryonic 
wind driven by neutrino heating above the neutrinosphere. 

In our multi-dimensional simulations the shock stagnates 
only for -120-130 ms and then expands to (Rghock) ~ 
6000 km at 400 ms postbounce (p.b.) compared to only 
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Fig. 2.— Mass-shell (radii for chosen enclosed masses) evolution of the 3D 
explosion with colors showing rms values of the nonradial velocity. The red 
solid line marks the angle-averaged shock radius, the dashed lines separate 
regions with mass spacings of 0.01, 0.1, 0.01, and 0.001 M© (from inside 
outwards). The 3D simulation was started from a ID model at ~ 10 ms after 
bounce (vertical dotted line). 


~2000km in the ID case (Fig.|^. Angle-averaged quantities 
X(r) are computed according to 


<A(r)) = 


f X(r)dQ 


(3) 


Positive total energies in the postshock layer develop only 
shortly after the onset of outward shock acceleration. Inter¬ 
estingly, in 3D the shock expands faster and the explosion 
energy increases more steeply and to a higher value than in 
2D. At fpb = 400 ms the diagnostic energy reaches 


£^fpi(400ms) ~ (0.77... 1.05) x 10™erg (4) 


in 3D with an increase at a rate of ~(0.75... 1.25)x 10™ erg s“* 
due to the PNS wind for the cases of minimal and maximal en¬ 
ergy, respectively, whereas it is only ~(0.67 ... 0.92) x 10^^ erg 
in 2D with a similar growth rate by the PNS-wind power as in 

3E0 

While the onset of the explosion in multi-dimensions is 
aided by buoyant convection, the steeper rise of the explo¬ 
sion energy in 3D can be understood by a smaller net energy- 
loss rate, Qcooh in the neutrino-cooling layer during 100 ms < 
fpb < 200 ms (Fig. 1^. Qcooi is evaluated between the mean 
gain radius, Rgain, and an inner radius Rq = = 3), 

which is (somewhat arbitrarily) defined at an optical depth of 
3 for Vg (all quantities are evaluated on angle-averaged pro¬ 
files). In contrast, the integrated net neutrino-heating rate 
in the gain layer, Qgain, and the heating efficiency, = 
2gain(^y, + (with Ey. being the direction-averaged lu¬ 

minosities of Ve and v^), do not exhibit any significant 3D-2D 
differences (Fig.|^. 

The reduced energy-loss rate in the cooling layer for 
100 ms < fpb < 200 ms is associated with a systematically 

lower specific internal energy, ^dVp^ j 

averaged over the cooling-layer mass Mcooi = dVp 

(Fig. [^. This suggests lower temperatures in the cooling 
layer, less efficient neutrino emission, and therefore an inward 


^ Tests with doubled angular resolution and different seed perturbations 
revealed ~4% variations of £^^^^(400 ms). 


shift of the gain radius in 3D compared to 2D: 

(Fig. [ 3 ]). 

Consequently, the mass of the gain layer, Mgain = 
r</?shock> jg larger and, correspondingly, the mass in the 

‘^^gain 

cooling layer, Mcooi, smaller in 3D. While the 3D-2D differ¬ 
ence of Mcooi varies with time because the definition of Rq 
depends on model-specific details, the gain-mass difference 
grows to AMgain ~ 1.2x10“^ Mq between 100 ms and 200 ms 
and remains essentially constant afterwards. The nuclear re¬ 
combination energy of this mass difference, 

AMgain X erecomb ~ (1-6... 2.0) X lo"'’ erg > (5) 

for ^recomb ~ (7 ... 8.8) MeV/nucleon (6.76... 8.49) x 
10^^ ergg“^ (depending on whether free nucleons recombine 
to Qf-particles or iron-group elements) can easily account for 
the 3D-2D differences of the explosion energy in Fig.[^ 
Turbulent pressure in the gain layer is understood t o oe sup- 
portive for the onset of neutrino-driven explosions ([Murphy 
jet al.|2013t[Couch & Qtt|20T^[Muller & Janka|2014bD . Ho^ 
ever, for our low-mass SN progenitor this effect cannot be 
crucial for the observed 2D-3D difference. Despite the faster 
and more energetic explosion, the 3D model exhibits a lower 
turbulent kinetic energy in the gain layer during the crucial 
period 100 ms < ^pb < 200 ms, visible by the nonradial con- 

tribution ^ 

into account the larger mass in the gain layer, it means that 
the nonradial fluid motions are much less vigorous in the 3D 
case. 

In contrast, the nonradial kinetic energy in the cooling layer, 
+ ^ 0 )’ smaller in the 2D simulation 
despite the higher mass of this layer in 2D. This suggests that 
nonradial motions are much more efficiently damped by dis¬ 
sipation in stronger deceleration shocks and due to symmetry 
constraints (reflective axis boundaries) in the cooling layer of 
the 2D model. 


3.2. Detailed Analysis of2D-3D Flow Differences 


The reason for more favorable explosion conditions in 3D 
becomes clearer through radial profiles of the turbulent accre¬ 
tion flow in the gain and cooling layers (Figs, ffl and |^. Here, 
angular averages of quantities X(r) over infallmg matter, i.e., 
over regions with radial velocity Vr < 0, are computed accord¬ 
ing to 




f X(r)&(-Vr)dQ 

f 0(-Vr) dn 


( 6 ) 


where 0(v) is the heaviside function with 0(v) = 1 for x > 0 
and 0(x) = 0 otherwise. 

Radial velocity profiles of angle-averaged over con¬ 

vective downdrafts in the postshock accretion layer, ex¬ 
hibit stronger local variations with higher extrema and pro¬ 
nounced intermittency-like behavior in 2D, whereas the 
angle-averaged infall velocities in the 3D model appear 
smoother and show less extreme fluctuations (Fig. left). 
This points to a larger number of smaller-scale convective 
structures, thus reducing the statistical variance, and suggests 
less vigorous radial mass motions in the 3D model, which is 
compatible with the lower nonradial kinetic energy in the gain 
layer between 100 ms and 200 ms p.b. (cf. Fig.[^. 

This conclusion is supported by cross-sectional cuts 
(Fig.§, where narrow convective downflows in 2D possess 
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Fig. 3.— Explosion parameters as functions of post-bounce time. ID results are displayed by dotted lines, 2D by thin solid lines, and 3D by thick solid lines. 
Angle-averaged shock radii (upper left)-, explosion energies with upper limits in red, lower limits in blue (see text) and 3D-2D differences with dashed lines 
(bottom lefty, kinetic energies of non-radial mass motions in the gain (red) and cooling layers (blue; upper middle)', average specific internal energies in the 
cooling layer (lower middle)', total net neutrino-heating rates in the gain (red) and cooling layers (blue) and neutrino-heating efficiencies (dashed, scaled by a 
factor of 10; top right)', masses in the heating (red) and cooling layers (blue) and difference of angle-averaged gain radii, -A/?gain = “(^gSn “ (dashed; 
bottom right). 


higher velocities, in agreement with radial profiles of av¬ 
eraged over inflows and time-averaged over 100 ms < fpb < 
200 ms (Fig.ffl. The higher infall velocities in 2D in the gain 
layer and a mr part of the cooling layer are associated with 
a greater mass-accretion rate around the gain radius (Fig.[^. 
Moreover, in 2D a larger fraction of the sphere in the gam 
layer is subtended by downfiows, which is expressed by the 
“downfiow filling factor”. 


f 0(-Vr) dQ 


(7) 


in Fig.pl In contrast, in the cooling layer 
these findings are in line with higher kinetic energies of ra¬ 
dial and nonradial flows in the gain layer and slightly lower 
nonradial kinetic energy in the cooling layer for the 2D model 

(FigM 

At first glance, the higher mass-accretion rate around /?gain 
in 2D seems to contradict the 2D-3D differences of the flow 
geometry visible in Fig. In 2D one counts one or two 
prominent convective downdrafts, consistent with a convec¬ 
tive cell pattern characterized by low spherical-harmonics 
modes, i.e., large-scale structures with few downfiows. In 
contrast, many more downflow funnels exist in 3D, although 
with smaller infall velocities. Nevertheless, the 2D model 
channels more mass per time through the gain radius. Ob¬ 
viously, the total volume of toroidal “sheets” formed by axi- 


symmetric (2D) downfiows is bigger than the volume encom¬ 
passed by 3D accretion funnels. Exactly this is expressed by 
the larger in the gain layer, which, together with bigger 
infall velocities, accounts for the higher mass-accretion rate 
in 2D. 


3.3. JD Turbulence Facilitating Explosion 

The reduced rate and less powerful infall of mass through 
the gain radius is therefore the fundamental reason for more 
favorable explosion conditions in 3D. Higher downflow ve¬ 
locities in 2D lead to more violent deceleration of the accre¬ 
tion downdrafts in shocks as they penetrate into the cooling 
layer (see Fig. ^and the more dramatic decrease of | {Vr)y^<o \ 
in 2D in Fig.|4 ). On the one hand this causes stronger gravity- 
wave activit}^elow which can be recognized by larger 
angular variations of Vr for the 2D model in this region in 
Fig.[^ On the other hand the impact of downfiows dissipates 
kinetic energy, for which reason the 2D simulation shows a 
lower kinetic energy of nonradial mass motions in the cooling 
layer despite its higher nonradial kinetic energy in the heat¬ 
ing layer (Fig.[^. More important, however, is the enhanced 
shock heating m the cooling layer, which leads to higher tem¬ 
peratures in this region in 2D (Fig.|^. Estimates confirm that 
this difference of the angle-averaged temperature of several 
10^ K (or several percent of the local angle-averaged temper¬ 
ature) is responsible for the larger net cooling rate of the 2D 
model just below /?gain (Eig. because local temperature 
fluctuations are considerably bigger than the relative differ- 
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Fig. 4.— Profiles of angle-averaged radial velocities of infalling (v^ < 0) matter, outgoing shock (sharp blue-white discontinuity), mean gain radius, Rgain, and 
inner boundary of cooling layer, Rq, as functions of post-bounce time for 2D {top left) and 3D models {bottom left). Time-averaged (over 0.1 s< fpb < 0.2 s) 
radial profiles of angle-averaged downflow velocities, (vr>v^<o» of matter with Vr < 0; upper middle)', mass-infall rates in downflows, Mv^<o {bottom middle)', 
“surface filling factors” of downflows, crv^<o {upper right)', and angle-averaged specific net neutrino-heating/cooling rates, {q), and 3D-2D temperature difference, 
A {T) = <r>^^ - {bottom right). Thin solid lines show 2D, thick solid lines 3D results, the vertical dotted lines indicate the time-averaged gain radius in 
3D. 



Vr {r > -Rgain) [ 10 ® CHI S Vr {v < -Rgain) [ 10 ® CHI S 

Fig. 5.— Cross-sectional cuts with radial velocities in the x-z plane at 0.16 s {left) and 0.18 s after bounce {right) for 2D {left semi-panels) and 3D models (right 
semi-panels). The mean gain radius, /?gain, is marked by a dashed circle. For better visibility different color scales are used outside and inside of /?gain (left and 
right color bar, respectively). 
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ences of the angular averages and neutrino emission rates by 
^--captures and ^^^“-pair annihilation scale with and 


respectively. Consequently, one finds during 

170 ms < ^pb < 220 ms (Fig.[^ bottom left). Since the temper¬ 
ature in the gain layer is only marginally higher in 3D, the net 
heating rates between 7?gain and shock are hardly distinguish¬ 
able for 2D and 3D (Fig.0 bottom right). 

The 2D-3D differences m the convective postshock layer, in 
particular higher downflow velocities and larger mass-infall 
rate in 2D, can be attributed to well-known differences of 2D 


and 3D turbulence discussed previously (e.g., 

Hanke et al. 

20 12[ [Burrows et al.|2012[ 

Dolence et al.|2013 

Couch|2013l 

Couch & O’Connor|2014l[] 

Dolence et al.|2013[ 

Abdikamalov[ 

et al. 
tation 

2014J: Turbulent energy cascading leads to fragmen- 
of large-scale flows to smaller-scale vortices in 3D, 


whereas the inverse energy cascade in 2D tends to enhance 
kinetic energy on the largest possible scales. The turbulent 
fragmentation (e.g. through Kelvin-Helmholtz instability) of 
the postshock accretion flow affects the flow characteristics 
most strongly near the gain radius, for which reason 2D-3D 
differences of Mv,<o increase as the infalling gas approaches 
^gain (Fig. 0- Because of enhanced turbulent energy in small- 
scale vortices, the 3D flow becomes less coherent, and tur¬ 
bulent small-scale structures keep mass in the gain layer and 
reduce Mvr<o^ while in 2D the mass-infall rate continues to 
grow all the way towards 7?gain. 


4. CONCLUSIONS 

Performing simulations of neutrino-driven explosions for 
a 9.6 Mq star with the Prometheus-Vertex code in ID and 
multi-D, we found significant 2D-3D differences of the accre¬ 
tion flow in the convective postshock layer. Turbulent frag¬ 
mentation to small-scale vortices keeps more matter in the 
gain layer in 3D and decreases the rate and velocity of mass 
accretion into the cooling layer. This reduces the dissipation 
of kinetic energy in the cooling region, lowers the temper¬ 
ature there and decreases the energy loss by neutrino emis¬ 
sion. The corresponding inward shift of the gain radius leads 
to an increase of the mass in the gain layer, whose higher nu¬ 
cleon recombination energy accelerates shock expansion in 
3D and boosts the explosion energy by >10% to >10^^ erg 


with a steeper rise at the beginning of the explosion. 

We thus presented the first successful neutrino-driven ex¬ 
plosion in a fully self-consistent 3D simulation and, for 
the first time, identified effects that foster and strengthen 
neutrino-driven explosions in 3D compared to 2D. In contrast 
to the unfavorable influence of 3D turbule nce observed in pre- 
vious studies of more massive stars (e. g., |Hanke et al. 2012 


Couch|20T3| [Couch & Q’Connor|2014tlTakiwaki et al.]20l'^ 


2014[|Abdikamalov et al.|2014j ), the decrease of the mass in 

flow rate through the gain radius and reduced neutrino-energy 
loss in the cooling layer below 7?gain have a healthy influence 
on the explosion in our 9.6 Mq progenitor in 3D. 

It is difficult to speculate whether these effects could play 
an important role also for more massive stars, in parti cular 
when convection dominates the post shock dynamics ( |Bur-| 
rows et aL||2012[ Murphy et al.||2013| ) and curr ent numerical 


resolutioiTdeficiencies (see, e^, discussions byjRadice et al. 


|2Q15t|Abdikamalov et al.|2Q14t|Hanke et al.|2Q12j ) have been 

overcome. We emphasize that our present grid resolution is 
insufficient to achieve convergence and to capture fully de¬ 
veloped turbulence. However, we attribute our 2D-3D differ¬ 
ences to generic differences of convective downflow dynam¬ 
ics, which are likely to remain valid even for higher resolu¬ 
tion. We have identified turbulent fragmentation of 3D flows 
in the gain layer as helpful for stronger explosions than in 2D, 
in sharp contrast to disadvantageous consequences of 3D ef¬ 
fects for the onset of explosions reported for more massive 
progenitors in previous works. Higher resolution might even 
enhance our described effects and their explosion-facilitating 
consequences rather than endangering them. 
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